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Abstract 

Many front-end applications of electron linear accelerators rely on the production of temporally-compressed bunches. 
The shortening of electron bunches is often realized with magnetic bunch compressors located in high-energy sections 
of accelerators. Magnetic compression is subject to collective effects including space charge and self interaction 
via coherent synchrotron radiation. In this paper we explore the application of magnetic compression to low-energy 
(~ 40 MeV), high-charge (nC) electron bunches with low normalized transverse emittances (< 5 fim). 

Keywords: Photoinjector, Linear accelerator, Electron beam. Magnetic bunch compressor. Space charge. Coherent 
synchrotron radiation. Flat beam 



1. Introduction 

Most of the photoinjectors being used for generation of 
bright electron bunches for, e.g., free-electron laser (EEL) 
applications consist of generating and rapidly accelerating 
the electron bunch to high energy and subsequently short- 
ening the bunch using a magnetic bunch compressor fH]. 
The only deviation to such a design is the combined ac- 
celeration and compression using velocity bunching |2]. 
Attempts to operate low-energy magnetic bunch compres- 
sors have to date been inconclusive |3] or deemed incom- 
patible with the production of low-emittance beams I'Tj. 
In this paper we explore and demonstrate via numeri- 
cal simulations that low-energy bunch compression per- 
formed on a ~ 40-MeV electron beam can be viable de- 
pending on requirements. We especially present trade- 
off curves between transverse emittance and peak current 
for several cases of electron-bunch charge. In addition, 
our simulations are performed with several computer pro- 
grams thereby enabling a benchmarking of very different 
approaches for modeling collective effects and especially 



coherent synchrotron radiation (CSR) ||5l |6|]. Our study 
considers the magnetic bunch compression planned in the 
40-50 MeV photoinjector of the Advanced Superconduct- 
ing Test Accelerator (ASTA) currently under construction 
atEemilab [7,1]. 



2. Accelerator beamline overview 

The compression of a ~ 40-MeV electron bunch via 
magnetic compression is investigated for the case of the 
ASTA photoinjector diagrammed in Eig. [T] The beam- 
line includes a photoemission electron source consisting 
of a cesium telluride (Cs2Te) photocathode located on 
the back plate of a Ih-1/2 cell radiofrequency (RE) cav- 
ity operating at 1.3 GHz \Qi- The cathode is illuminated 
with a 3-ps ultraviolet laser pulse with uniform radial dis- 
tribution and a Gaussian temporal profile. The RE gun 
is surrounded by two solenoidal lenses that control the 
beam's transverse size and emittance. Downstream of 
the RE gun, the typical beam energy is ~ 5 MeV. The 
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bunches are further accelerated up to 50 MeV by two 
1.3-GHz superconducting RF (SCRF) accelerating cav- 
ities (labeled as CAVl and CAV2 in Fig. [B- A third 
SCRF cavity (CAV39) operating at 3.9 GHz will even- 
tually be incorporated to correct for nonlinear longitudi- 
nal phase space distortions |l3- 12]. Because of its su- 
perconducting nature, the ASTA facility produces elec- 
tron bunches repeated at 3 MHz arranged in a 1-ms 5- 
Hz RF macropulse. The downstream beamline includes 
quadrupoles, steering dipole magnets, and diagnostics sta- 
tions. A skew-quadrupole channel can be set up as a 
round-to-flat-beam transformer (RFBT) to convert an in- 
coming angular-momentum-dominated beam into a flat 



beam with high transverse emittance ratio [ 13l |14|] . The 
beamline also incorporates a four-bend magnetic bunch 
compressor (BCl) which, consists of four 0.2-m rectan- 
gular dipoles (Bl, B2, B3, B4) with respective bending 
angles of {+,-,-,+) 18°. The longitudinal dispersion of 
BCl is 7?56 = -0.19 m. Finally a single-shot longi- 
tudinal phase space diagnostics combining a transverse- 
deflecting cavity (TDC) with a vertical spectrometer will 
be instafled 1151 
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inal betatron functions downstream of CAV39 up to the 
cryomodule entrance is plotted in Fig.|2] 




Figure 2: Evolution of tlie liorizontal (blue) and vertical (red) betatron 
functions through the ASTA injector. The green rectangles indicate the 
location of quadmpole and dipole (smaller rectangles) magnets. The 
BCl compressor is located at s e [11.9,15.1] m. The origin of the 
horizontal axis (s = m, not shown) corresponds to the photocathode 
surface, (color online) 



22.0 3. Modeling methodology 



Figure 1: Injector configuration at ASTA. The "RF gun", "LI" and "L2" 
respectively correspond to the gun cavity and surrounding solenoid mag- 
nets, "CAVl", "CAV2", and "CAV39" are superconducting RF cavities. 
"RFBT" is the round-to-flat beam transfomier, and "BCl" refers to the 
magnetic bunch compressor, and Bl-4 are the dipoles of the chicane, 
with distance between the dipoles marked in the figure. The number 
below the beamline indicates the axial positions in meters w.r.t. the pho- 
tocathode surface. 

The beam dynamics through CAV2 were simulated 
with ASTRA and optimized using a genetic optimizer for 
several cases of charge and photocathode drive-laser con- 
figurations; see Ref. Iil6il . The resulting phase space 
distributions are used as a starting point for transport 
and compression through the beamline downstream of 
CAV39. The quadrupoles settings were optimized for 
the various operating charges using the single-particle dy- 
namics program elegant [17]. The evolution of the nom- 



The bunch-compression performance of BCl was ex- 
plored via numerical simulations using impact-z [18] and 
csRTRACK ifigll . Both programs model the beam as an 
ensemble of interacting macroparticles and integrate the 
equations of motion to advance the macroparticles along 
a user-specified beamline. 

In MPACT-z the space-charge (SC) interaction is mod- 
eled using a mean-field quasi-static particle-in-cell (PIC) 
algorithm and the (point-like) macroparticles are ad- 
vanced through the beamline using high-order transfer 
maps. Each beamline element is segmented into axial 
slices modeled by transfer maps. Between each transfer- 
map segment, impact-Z applies a space-charge "kick" 
evaluated from the mean-field PIC SC algorithm I.18l] . 
CSR eff'ects are included in mpact-Z using the one- 
dimensional formulation described in Ref. [20]. The one- 
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dimensional model is valid provided |121 



2)(s) <K 1, with = 
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where R(s) is the trajectory's radius of curvature and (Tx(s) 
and (Tzis) are respectively the transverse and longitudinal 
root-mean-square (RMS) sizes at the curvilinear beamline 
position s. 

In iMPACT-z the longitudinal charge distribution needed 
for the ID CSR model is obtained from a longitudinal bin- 
ning of the macroparticle ensemble. Convergence studies 
were carried out in order to determine the optimal num- 
ber of longitudinal bins, N^, to be used for both the SC 
and CSR calculations. Low values of A^, generally under- 
estimated the peak-current and therefore the collective ef- 
fects, while large values of N, introduce numerical noise 
that can lead to artifacts (e.g. numerically-induced mi- 
crobunching) l(22[l . The convergence study [23] revealed 
an appropriate value of Nz = 256 for a bunch represented 
by = 2 X 10^ macroparticles. The number of bins in 
the transverse dimensions was set to Nx - Ny = 16. The 
iMPACT-z simulations presented in the rest of this paper use 
this set of parameters. 

The program csrtrack was specifically developed to 
efficiently simulate the impact of bunch radiative self- 
interaction via CSR. Csrtrack incorporates several mod- 
els including a 2D particle-to-particle (P2P) model that 
directly computes the forces on macroparticles from the 
Lienard-Wiechart potentials evaluated at retarded times. 
These calculations are self-consistent and enable the com- 
putation of both the transverse and longitudinal force 
contributions from SC and CSR effects. Since the P2P 
model is computationally intensive, (the calculation time 
scales as A^^), csrtrack also includes an improved one- 
dimensional model referred to as the ID Projected (IDP) 
model. The IDP model uses the ID projection of the 
smoothed charge distribution convoluted with a kernel 
function ||24ll . Compared to the model of Ref. ll20ll . the 
IDP model is not limited to the ultra-relativistic regime. 
Csrtrack's P2P model treats each macroparticle as a 3D 
Gaussian charge distribution (referred to as "sub-bunch") 
in the (x, y, z) space with distribution 



g(x,y,z) = 



1 
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where cr/„ cr,, and cr|| are respectively the horizontal, ver- 
tical and longitudinal RMS sizes of the sub-bunches. The 
resulting beam's spatial charge distribution is <l>(x,y, z) = 
2^=1 Qjgix - Xj,y - yj,z - Zj) where Qj and (xj,yj,Zj) 
are respectively the i-th sub-bunch charge and LPS coor- 
dinates. In CSRTRACK, cr,, and o-y may be defined relative to 
the vertical and longitudinal RMS bunch sizes, cr,, and 0-,, 
respectively, and are adjusted along the bunch compres- 
sor as the dimensions of the bunch change. Due to the 
computational intensiveness of the P2P model, only 10* 
sub-bunches were used compared to 2 x 10^ used with the 
IDP model. This relatively-low number of sub-bunches 
requires a large o-|| of 5% of the longitudinal bunch length 
cr, for the IDP simulations and 10% for the P2P simula- 
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tions; for a detailed study see Ref. 112311 . The P2P model 
requires parameters for the sub-bunches in the horizontal 
and vertical dimensions, o-/, and cr,,. We chose cr,,= 0. lOcr,, 
and cr|p O.lOcr,. However, csrtrack does not allow cr/, to 
be set as a variable of the horizontal RMS width, so we in- 
stead chose 0.1 mm, which is on the order of 10% of the 
RMS size cr j. for all four of the bunch charges presented 
here. Both of csrtrack's models neglect collective forces 
in the vertical dimension. 

The ASTRA simulations of the beam generation and ac- 
celeration simulated up to 5 - 8.2 m from the photo- 
cathode were used as a starting point for our simulations. 
A PYTHON program, gluetrack 1,25,1 , was used to manip- 
ulate the beam distributions and generate macroparticle 
distributions suitable for impact-z and csrtrack. For all 
the bunch-compressor studies, impact-z was used to track 
the bunch distribution from s - 8.2 up to ,s = 11.8 m 
corresponding to 0. 1 m upstream of the entrance face of 
the first dipole (Bl) of BCl. This distribution was used 
as a starting point for the bunch-compressor simulations, 
which include BCl and a 1 .0-m downstream drift to allow 
for the influence of SC possible transient CSR effects. 

Simulations were performed for four cases of bunch 
charges ranging from 3.2 nC to 20 pC. For each charge 
the transverse emittance was optimized with astra ri6'l. 
The distributions were manipulated using gluetrack to 
adjust several parameters including their Courant-Snyder 
(C-S) parameters, and longitudinal phase space (LPS) 
chirp C = -{Zi6i)lcrl. where (z,, 6i) are the coordinates in 
the LPS, the (m) indicates the statistical averaging of vari- 



able u over the LPS distribution and cr^ ,■ = {z.) 
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also modeled the effect of CAV39 by numerically remov- 
ing the second order correlation in the LPS distribution. 



Table 1: Transverse and longitudinal beam parameters 0.1 -m upstream 
of Bl dipole entrance face. Only the Courant-Snyder parameters were 
fixed while the other parameters depend on the bunch charge or upstream 
beamline settings. 



Parameter 



Value Units 



0!x,i 

c 

total energy 



3 

1.6 

-1.6 
[1.0,6.0] 
38.6 



m"' 
MeV 



The charge-independent beam parameters computed 
0.1 -m upstream of dipole magnet Bl are summarized in 
Tab.[T] The beam distribution was matched to achieve the 
transverse C-S parameters listed in Tab. [T](see discussion 
below) upstream of BCl . The initial LPS chirp was tuned 
by removing the second-order correlation and scaling the 
first-order correlation between 5, and z,. The other LPS 
parameters and emittances (shown for the four charges in 
Tab. IJI are inherent to the generation process and were 
not adjusted. The initial LPS distribution for each of the 
four charges appears in Fig. |3]with its linear correlation 
removed. The S -shaped LPS is a remnant of space charge 
effects during the bunch generation and transport before 
acceleration in CAVl and CAV2 |26]. As expected larger 
charges yield higher total fractional momentum spread. 



Table 2: Initial normalized transverse £v/v,i ™d longitudinal £; , emit- 
tances and RMS bunch length o". , for the four cases of charge considered 
in this paper The parameters are computed 0.1 -m upstream of dipole 
magnet B I's entrance face. 
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Figure 3: LPS distributions 0. 1-m upstream of dipole magnet Bl's en- 
trance face for 3.2 (a), 1.0 (b), 0.25 (c) and 0.02 nC (d) bunches. The 
distributions were obtained from simulations of the photoinjector beam 
dynamics in astra; see Ref. 1 16]. The ordinates z > correspond to the 
head of the bunch. 



Simulations were performed for LPS chirps C e 
[1 .0, 6.0] m"' with the bulk of the simulations performed 
around the maximum-compression value C - -l/Rsd - 
5.2 m ', corresponding to enhanced collective effects. 

The initial values for the betatron functions were se- 
lected such that the beam experiences a waist between 
the third and fourth dipole i27ll . Simulations performed 
with the IDP csrtrack model also confirmed there is a set 
of horizontal C-S parameters that minimizes the bending- 
plane emittance growth as displayed in Fig. |4] The up- 
stream magnets were tuned to provide the incoming hor- 
izontal C-S parameter (/Sv./, a^j) - (8.0 m, 3.0) shown in 
Fig. |2] Lastly, the final RMS bunch lengths as a func- 
tion of the initial LPS chirp are given in Fig. |5] as ob- 
tained using iMPACT-z's SC+CSR model. In addition, there 
are slight variations in the bunch length between the four 
different models due to their collective effects inhibiting 
compression to varying degrees. 



4. Benchmarking of numerical models 

The beam dynamics simulations throughout BCl were 
performed for several degrees of bunch compression (con- 
trolled with the LPS chirp) for the four cases of bunch 
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Figure 4: Contour plot of the final nomialized horizontal emittance {£■, 
in fim) as a function of the C-S parameters /3i j and cv, , 0. 1-m upstream 
of dipole Bl. The simulations were performed with csrtrack's ID- 
Projected model and a bunch charge of 3.2 nC . From this data, we 
selected the values for our simulations, (ySv,,, as,i) = (8.0 m, 3.0). (color 
online) 
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Figure 5: RMS bunch length cr- downstream of BCl asafunction of the 
LPS chirp for various bunch charges using Impact-Z's CSR+SC model 
model, for 3.2-nC (blue), 1.0-nC (red), 250-pC (green), 20-pC (ma- 
genta) bunch charges. The inset plot corresponds to a close-up around 
chirp values that achieve minimum RMS bunch lengths, (color online) 



charges that will be used at the ASTA. Four simulation 
algorithms were used as described earlier Because of 
its computational effort, crtrack's P2P model simulations 
were restrained to a smaller range of values for C. 

The LPS 1.0-m downstream of the B4 dipole simu- 
lated with iMPACT-Z (SC -I- one-dimensional CSR mod- 
els) and csRTRACK (P2P model) are shown in Fig. |6] [(a)- 
(d)] and [(e)-(f)] respectively for the four charges listed 
in Tab. |2l Despite the vastly different algorithms used by 
these two programs, the LPS distributions displayed very 
similar distortions including those at the small-scale lev- 
els. FigureQsummarizes the evolution of peak current as 
a function of the initial LPS chirp for the four numerical 
models. Likewise the longitudinal emittances computed 
with iMPACT-z (SC-hCSR) and csrtrack (P2P) for the case 
of maximum compression are in decent agreement; see 
Tab.m 

The transverse-emittance after compression is shown 
in Fig. [8] Csrtrack's P2P model consistently offers 
the greatest emittance growth, followed by impact-z's 



Table 3: Final normalized longitudinal £, at maximum compression 
(C = 5.2 m"') simulated with impact-z (SC+CSR) and csrtrack (P2P). 





IMPACT- Z 


CSRTRACK 




SC-hCSR 


P2P 


Q(nC) 


e- (//m) 


B, (//m) 


3.2 


267 


261 


1.0 


118 


105 


0.250 


61.5 


57.8 


0.020 


10.5 


11.6 



SCh-CSR model, as these are the only two models that 
account for both SC and CSR effects. However, the 
P2P model also includes transverse CSR forces and a 
more elaborate model for longitudinal CSR. Our previ- 
ous study fl^] showed that the influences on final emit- 
tance from using too-few macroparticles as well from as 
the randomization in the down-sampling of initial dis- 
tributions were both much smaller than the discrepancy 
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Figure 6: LPS at BCl exit for impact-z's (a-d) and csrtrack's (e-li) 3D 
models, for 3.2-nC (a,e), 1-nC (b,f), 250-pC (c,g), and 20-pC (d,h) buncli 
cliarges, zoomed in to show details, for C = 5.2 m"'. (Red line) longi- 
tudinal current projection, with arbitrary scale and offset. Note that the 
horizontal and vertical axis ranges are different for each plot. The ordi- 
nates z> correspond to the head of the bunch, (color online) 



between the SC-i-CSR and P2P models. The emittance 
growth observed from csrtrack's ID and impact-z's SC- 
only model, indicates that CSR accounts for most of the 
emittance dilution at higher charge. For the low-charge 
simulations {Q = 250 and 20 pC) the relative importance 
is reversed, with SC contributing more to the emittance 
degradation than CSR. Of the models presented here, only 
impact-z's include SC in both the vertical and the horizon- 
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Figure 7: Peak currents / versus energy chirp for impact-z's SC+CSR 
(green), impact-z's SC (blue), csrtrack's ID? (red), and csrtrack's P2P 
(magenta) models, for 3.2-nC (a), 1.0-nC (b), 250-pC (c), and 20-pC (d) 
bunch charges, (color online) 



tal planes, and are shown in Fig.|9] As vertical emittance 
growth is entirely the result of SC, the inclusion of impact- 
z's CSR model reduces vertical emittance growth due to 
the reduced compression. 

Finally, the evolution of the slice parameters during 
compression were explored. For this analysis the beam 
is divided into axial slices of equal longitudinal length 
6z - 20 fim. A statistical analysis on the population con- 
tributing to each slice was performed to yield the slice 
emittances, energy spread and peak current. A com- 
parison of the slice bending-plane emittance and energy 
spread between impact-z and csrtrack's P2P model of 
slice-emittance as the chirp is varied appears in Fig. [TO] 
The level of agreement between the two programs is of 
the same order as what observed for the bunch parameters. 
Fig. nn summarizes the evolution of the slice horizontal 
emittance in the slice with the highest peak current within 
the bunch. When the beam is greatly over-compressed 
(C = 6.0 m '), the LPS may be double-peaked, and the 
transverse brightness as defined here may not be an appro- 
priate measure of the bunch's utility. The curve demon- 
strates how little slice emittance growth occurs for partial 
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Figure 8: Final liorizontal emittances for eacli of tlie different bunch 
cliarges with impact-z's SC-only model (blue), csrtrack's ID CSR- 
model (Red), impact-z's SC+CSR model (Green), and csrtrack's P2P 
model (magenta), for 3.2-nC (a), 1-nC (b), 250-pC (c), and 20-pC (d) 
bunch charges, (color online) 



Figure 9: Final vertical emittances for each of the different bunch 
charges with impact-z's SC-only model (blue) and impact-z's SC+CSR 
model (Green), for 3.2-nC (a), 1-nC (b), 250-pC (c), and 20-pC (d) 
bunch charges. Csrtrack does not compute vertical forces, so the emit- 
tance remains roughly constant along the bunch compressor, (color on- 
line) 



compression. 

5. Expected beam dilution and trade-offs 

A large number of accelerator applications re- 
quire beams with high-peak-currents and low-transverse- 
emittances. These requirements conflict with each other 
as collective efifects, which dilute the beam's phase space 
and emittances, increase with peak current. A commonly- 
used figure of merit is peak transverse brightness = 



^ 112811 . Figure [12] summarized the evolution of 
as function of the LPS chirp for the four cases of bunch 
charges. The figure combines the data provided in Fig. |7] 
and Fig. [8] Despite the lower-charge bunch result in 
smaller peak current at lower charges (see Fig. 0, the 
transverse brightness increases with lower bunch charges. 
The main factor at play in this reduction is the lower initial 
transverse emittances, and more importantly, the reduced 
dilution of the transverse emittances during compression 
in BCl due to the weaker collective effects (CSR and SC). 



In addition we note that the maximum achieved trans- 
verse brightness does not necessarily occur at maximum 
compression. This is due to the larger peak currents 
at maximum compression that drive collective effects, 
wherein the relative emittance growth driven by collective 
effects is greater than the relative increase in peak current, 
a trade-off similar to that of going down to lower bunch 
charges. A summary of the achieved maximum value of 
B^ appears in Fig. [13] 

The trade-off between obtained peak current and is 
shown in Fig. [141 only data associated to LPS chirp up 
to maximum compression at C < 5.2 m ' are displayed, 
as over-compression results in lower peak currents with 
generally larger emittance dilutions. 

The simulation data points to several conclusions about 
the parametric trade-offs that must be considered. First, 
the emittance growth, particularly the slice emittance, is 
greatly reduced at lower degrees of compression, partic- 
ularly C < 4 m \ which corresponds to compression 
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Figure 10: Example of final normalized transverse slice emittances (top 
row) and energy spread (bottom row) evolution within a 1-nC bunch for 
four cases of compression C = 4.0, 5.0, 5.25, and 6.0 m"' respectively 
shown as blue, green, red and turquoise traces). Plots (a) and (c) cor- 
respond to iMPACT-z simulations while plots (b) and (d) are results from 
csrtrack's P2P model. Emittance and energy spread values associated 
to slices that contain too-few number of macroparticle for meaningful 
statistical analysis are set to zero. The heads and tails of the bunches 
are sparsely populated (see Fig.|6]for reference), particularly for the P2P 
simulations which use only 5 % of the number of particles used in the 
iMPACT-z simulations. 



to around one-third of the initial bunch length. Second, 
using lower bunch-charges is preferred for experiments 
that require high transverse brightnesses, due to the lower 
emittance growth from collective effects justifying the 
lowered peak current. For the 20-pC bunch charge, the 
regime with C < 4.0 m ' results in horizontal emittance 
growth that is under 10% of the initial horizontal emit- 
tance, regardless of which of the simulation codes is used. 
Under-compression is discussed further in section 6. 1 . 

6. Applications 

The main motivation that led to the inclusion of BCl in 
the ASTA's photoinjector is to provide a weak compres- 
sion necessary to avoid significant energy spread to be ac- 
cumulated during acceleration in subsequent cryomodule. 



Figure 1 1 : Final normalized transverse shce emittances « in the slice 
with the highest peak current computed with impact-z's SC-only model 
(blue), csrtrack's ID CSR-model (red), and impact-z's SC+CSR model 
(green), and csrtrack's P2P model (magenta), for 3.2-nC (a), 1-nC (b), 
250-pC (c), and 20-pC (d) bunch charges. 



It is anticipated that a second-stage bunch compressor will 
eventually be installed at higher energy to compress the 
bunch to high peak current while mitigating phase-space 
dilution. However, in light of the studies presented in the 
previous Sections it is worth investigating other possible 
applications of the BCl compressor as discussed below. 

6.1. Multi-stage bunch compression 

One of the motivations for exploring low-energy bunch 
compression is to moderately compress the bunches be- 
fore injection in an accelerating structure and then further 
compress at high energy [4]. In Ref. [4], the low-energy 
bunch compression was accomplished at ~ 15 MeV and 
resulted in intolerable beam degradation at the Tesla Test 
Facility I (TTFl) and consequently dismantled from the 
beamline during its upgrade as the Free electron LASer in 
Hamburg (FLASH) user facility. The requirement on the 
first-stage bunch compression is to provide bunch lengths 
that satisfy cr- <k A/(2n) (where A is the wavelength of 
the RF wave associated with the subsequent accelerating 
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LPS Chirp (m" LPS Chirp (m" 

Figure 12: Peak transverse brightness ^± = 5^27^ versus energy chirp 
for iMPACT-z's SC+CSR (green), impact-z's SC (blue), csrtrack's IDP 
(red), and csrtrack's P2P (magenta) models, for 3.2-nC (a), 1.0-nC (b), 
250-pC (c), and 20-pC (d) bunch charges, (color online) 
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Figure 13: Maximum peak transverse brightness fix = 4^2^ ^ versus 
bunch charge for impact-z's SC+CSR (green), impact-z's SC (blue), csr- 
track's IDP (red), and csrtrack's P2P (magenta) models. Each data 
point is a maximum from each line in Fig. 1121 (color online) 



Structure). The latter condition insures that no significant 
LPS quadratic distortion is imparted during acceleration 
in the subsequent linac. For the L-band linac of ASTA 
this sets the upper requirement cr- < 800 fim |4]. For 
a Gaussian distribution this requirement corresponds to a 
/ ^ 500 A at g = 3.2 nC and C ~ 3.7 m ' resulting in 
the tolerable bending-plane emittance dilution as shown 
in Tab.|4]along with compression to cr^ - 800 jum for each 
of the other bunch charges. Relative emittance growth is 
significantly smaller for 1-nC bunch charge and below. 

6.2. High peak current production 

Despite their relatively poor transverse emittance, fully 
compressed bunches could be used to generate copious 
amounts of radiation via a given electromagnetic process. 
The spectral-angular fluence emitted by a bunch of » 1 
electrons from any electromagnetic process is related to 
the single-electron spectral fluence, via 

L,^ \AN + N^\S{(jj)\\ (3) 



where o) = 2nf (f is the frequency) and S (w), the bunch 
form factor (BFF), is the intensity-normalized Fourier 
transform of the normalized charge distribution S(t) I2M- 
The former equation assumes the bunch can be approxi- 
mated as a line charge distribution and is practically valid 
as long as the rms bunch duration cr, and transverse size 
satisfy ctj. <k ccrt/y where y is the Lorentz factor and 
c is the velocity of light. When the BFF approaches unity, 
4-^1^, oc A^^ and the radiation is termed "coherent radia- 
tion". 

At ASTA the availability of a superconducing linac 
coupled with a non-interceptive radiation-generation 
mechanism (e.g. diffraction radiation [30]) could lead to 
the production of single-cycle THz pulses repeated at 3 
MHz over 1-ms. As an example we consider the worst- 
case scenario of a fully compressed 3.2-nC bunch; the de- 
pendency of the BFF over frequency appears in Fig. [15] 
(left plot). The BFF starts to take off" at frequency lower 
than / - 1 THz thereby supporting the generation of 
coherently-enhanced radiation at these frequencies. A 
limitation might come from the large transverse emittance 
that would prevent the beam to be focused to a transverse 
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Figure 14: Final normalized transverse emittances £j versus peak cur- 
rents / for LPS chirps using impact-z's SC+CSR (green), impact-z's SC 
(blue), csrtrack's IDP (red), and csrtrack's P2P (magenta) models, for 
3.2-nC (a), 1.0-nC (b), 250-pC (c), and 20-pC (d) bunch charges. Only 
data corresponding to chirp values C e [1.0,5.2] m"' are displayed, 
(color online) 



spot RMS size below the required cr^ly value. However 
a statistical analysis indicates that the central part of the 
beam containing approximately 15% of the beam popu- 
lation (or 500-pC out of the original 3.2-nC bunch) has 
emittances below 10 ^va resulting in beam cr^ < 100 ^va; 
see Fig. [15] (right plot). It should be pointed out that the 
lower-charge cases investigated in the previous section 
would result in shorter pulses with associated BFFs that 
contain higher-frequency content (see also Fig.|5). 

6.3. Compressed flat-beam generation 

An important asset of the ASTA photoinjector is its 
capability to generate beams with high-transverse emit- 
tance ratios known as flat beams. Immersing the photo- 
cathode in a magnetic field introduces a canonical angu- 
lar momentum (L) - eBocrl, with Bq the magnetic field 
on the photocathode surface, and cTc the RMS transverse 
size of the drive-laser spot on the photocathode |31]. As 
the beam exits the solenoidal field provided by lenses LI 
and L2, the angular momentum is purely kinetic result- 
ing in a beam coupled in the two transverse planes. Three 



Table 4: Final normalized transverse and longitudinal £; emittances 
for the four cases of charge considered in this paper. The initial LPS 
chirp (C) was optimized for each charge to yield a final bunch length 
CTj ^ 800 /im, based on the scan presented in Fig. [5] The simulations 
were performed with impact-z's SC+CSR model. The values displayed 
in this Table should be compared with the pre-compression values sum- 
marized in Tab. |2] 



Q(nC) 


C(m-') 


e V (yum) 


Ey (jjm) 




3.2 


3.7 


10.56 


4.24 


88.7 


1.0 


3.2 


3.03 


1.98 


34.2 


0.250 


3.1 


0.623 


0.594 


14.7 


0.020 


1.9 


0.296 


0.293 


1.83 
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Figure 15: Bunch form factor associated to the 3.2-nC fully-compressed 
electron bunch (left plot) and final horizontal (blue), vertical (red), and 
longitudinal (green) normalized emittances (left vertical axis) and peak 
current (dashed black line, right axis) of the bunch within a selected 
transverse radius (right plot). These simulations were performed near 
maximum compression with C = 5.2 m"' and 3.2-nC. (color online) 



skew quadrupoles in the beamline can apply the torque 
necessary to cancel the angular momentum ||32[ l33ll . As 
a result, the final beam's transverse emittance partition is 
given by 

/ £.2 



(4) 



where e„ is the normalized uncorrelated emittance of the 
magnetized beam prior to the transformer, /? and y the 
Lorentz factors, X = {L)l2p^, and p~ is the longitudi- 
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nal momentum. Note that the product e^.i^j,; = 

If compressed these flat beams may have appUcations in 

Smith-Purcefl FELs |34] or for beam-driven acceleration 

techniques using asymmetric structures I35ll . It may also 

be possible to mitigate the emittance growth in BCl by 

having a beam that is wide in the direction of the chicane 

bend. 

In this section, we explore the behavior of flat beams 
in the low-energy bunch compressor at ASTA, for the dif- 
ferent initial emittance ratios p s ex,i/ev,i- In order to 
produce these bunches, we took the 3.2-nC bunch pre- 
sented earlier and numerically scaled the macroparticle 
coordinates to produce the desired transverse emittance 
ratios while constraining the product Sx,iEy,i - 5^ jivc?. 
Due to the large transverse aspect ratio of the bunches, 
the criterion given in Eq{T]is generally not satisfied and 
it is therefore anticipated that the projected CSR model is 
inadequate, thus we use CSRtrack's P2P model to sim- 
ulate the flat beams and neglect the IDP model. The 
parameters used for flat beam simulations follow those 
used in the previous section, with the exception of the 
macroparticle horizontal size used in the csrtrack P2P 
model. Due to the much greater transverse dimension we 
set (Th = 0.2 mm. In addition, impact-z SC-hCSR sim- 
ulations were also performed to evaluate the emittance 
growth in the vertical plane. The simulated emittance 
growth is shown in Fig. [16] for a 3.2-nC bunch with an 
initial LPS chirp of C - 5.2 m As expected the rela- 
tive emittance dilution is reduced as the initial emittance 
ratio p increases. The agreement between csrtrack and 
iMPACT-z for the bending-plane emittance dilution is re- 
markable (within ~ 30 %) given the large transverse hori- 
zontal beam sizes. Impact-z predicts that the vertical emit- 
tance increases by a factor of 1.5 to 1.8 over the range of 
considered initial emittance ratios p € [1, 500]. The four- 
dimensional transverse emittance growth S4 = ye^ev is 
mitigated for the larger initial flat-beam emittance ratios. 

6.4. Double-bunch generation 

The production of shaped electron bunches has a large 
number of applications including the investigation of 
Wakefield and beam-driven acceleration techniques. Op- 
erating the low energy bunch compressor with LPS chirp 
C > 5.5 m"' leads to over compression and results in 
a structured longitudinal charge distribution. Figure [TTI 
confirms, for the case of 2 = 3.2 nC, that a bi-modal 




Figure 16: Bending plane transverse emittance ^ growth in BCl (red) 
simulated with csrtrack (dashed line) and impact-z (solid lines) as func- 
tions of the initial emittance ratio p = CoiTesponding impact- 

z results for the vertical emittance (^, blue solid line), and four- 
dimensional transverse emittance (^, magenta solid line), (color on- 
line) 



distribution could be generated with a separation between 
its peaks (~ 300 /vm) consistent with requirements from 
beam-driven acceleration such as plasma-wakefield and 
dielectric-wakefield acceleration techniques. In addition 
the distance between the peaks could be controlled to 
some degree by slight changes over the initial LPS chirp. 
The full-bunch and slice-at-peak-current horizontal emit- 
tance at C ~ 5.5 m"' are 67.0 and 75.1 fim, respec- 
tively, compared to 106 and 107 fim for the maximum- 
compression case (C ~ 5.2 m"'). These bending-plane 
normalized emittances of ~ 75 fim can still be focused to 
a sub-mm or sub-100-/im transverse spot size at respec- 
tively ~ 40 MeV and ~ 250 MeV (the latter energy cor- 
responds to acceleration of the 40 MeV beam into one of 
ASTA accelerating cryomodules). 
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z (mm) z(mm) 



lation and commissioning of the BCl bunch compressor 
along with the available diagnostics at ASTA will pro- 
vide a unique experimental platform for benchmarking 
the simulation codes currently available in regimes where 
the CSR and SC effects can play similar roles. 
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Figure 17: LPS distribution (gray colomiap in left plot) and current pro- 
jection (red trace and right plot) associated to over-compressed bunches 
with an incoming LPS chirp of C = 5.5 m"' . The ordinates z > corre- 
spond to the head of the bunch, (color online) 

7. Summary 

In this paper we presented numerical studies of a low- 
energy magnetic bunch compressor similar to the one be- 
ing installed at the ASTA facility at Fermilab. Our re- 
sults indicate that low-energy compression can be a vi- 
able path as a first stage compressor for a multistage com- 
pression scheme. In addition the capability of such a 
low energy compressor to provide high-peak currents at 
low energy (~ 40 MeV) could have important applica- 
tions when combined with the long-macropulse capabil- 
ity of the ASTA's superconducting accelerating module 
such as, e.g. the production of single-cycle THz radi- 
ation from diffraction radiation. As part of our investi- 
gation we used several computer programs and observed 
an acceptable agreement when simulating the evolution 
of the LPS distributions during compression. The sim- 
ulated transverse-phase-space parameters downstream of 
the bunch compressor have discrepancies that are inherent 
to the capabilities of each of the model but provide simi- 
lar emittance growth and trade-off curves. Based on these 
observations, simple (and faster) models could be used to 
optimize the bunch compression design with a quick turn 
around while first-principle model could provide high- 
fidelity simulations of the optimized designs. The instal- 
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